Molecular imaging using high-order harmonic generation and above-threshold 

ionization 
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Accurate molecular imaging via high-order harmonic generation relies on comparing the harmonic 
emission from a molecule and an adequate reference system. However, an ideal reference atom with 
the same ionization properties as the molecule does not always exist. We show that for suitably de- 
signed, very short laser pulses, a one-to-one mapping between high-order harmonic frequencies and 
electron momenta in above-threshold ionization exists. Comparing molecular and atomic momen- 
tum distributions then provides the electron return amplitude in the molecule for every harmonic 
frequency. We show that the method retrieves the molecular recombination transition moments 
highly accurately, even with suboptimal reference atoms. 

PACS numbers: 33.80.Rv, 42.65.Ky 
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When atoms or molecules arc irradiated by a strong 
laser field, high-order harmonic generation (HHG) takes 
place and high-frequency photons are emitted The 
interest in HHG from molecules is growing since observ- 
ing the radiation is a tool to investigate the structure 
of molecules [H-Q. The sensitivity of the emission spec- 
tra to the target structure can be understood within the 
three-step model, which provides a semiclassical interpre- 
tation of HHG in terms of (i) ionization, (ii) free prop- 
agation of the electron in the laser field and return to 
the parent ion, and (hi) recombination In good ap- 
proximation, the HHG intensity is proportional to the 
modulus squared of the recombination transition dipolc 
moment, or equivalently to the recombination cross sec- 
tion 0, 0; HI ■ When the electron continuum states are ad- 
ditionally approximated as plane waves, one can obtain 
molecular orbitals via a tomographic retrieval based on 
the Fourier transform of the HHG amplitudes measured 
from aligned molecules 0, HI • These HHG-based molec- 
ular imaging methods rely on the comparison of the har- 
monic emission to the one from a reference system with 
known electronic structure, typically a reference atom. 
Assuming that the molecule and the reference atom have 
the same properties concerning ionization probabilities 
and electron propagation, the recombination cross sec- 
tion of the molecule can be isolated. Clearly, no refer- 
ence atom with exactly the same ionization properties as 
the molecule exists. Therefore, a systematic method to 
correct for these deviations is highly desirable. 

Since the probability for recombination of a return- 
ing electron is very small, it is likely that the system 
remains ionized, and the electron can be detected as an 
above-threshold ionization (ATI) electron. ATI momen- 
tum distributions have also be used for molecular imaging 
It appears plausible to combine HHG and ATI to 
improve laser-based molecular imaging. For reasons out- 
lined in the following, however, no concrete method has 



been proposed up to now. 

Usually, multi-cycle laser pulses have been used to 
drive HHG and ATI. This means that many diffcrcrcnt 
electron trajectories can potentially contribute to the 
same harmonic frequency or to the same electron momen- 
tum. In the case of HHG, the same frequency is generated 
twice per optical half cycle, namely by the well known 
short and long trajectory [Tl(. In the case of ATI, the in- 
terference of contributions from two ionization times has 
been termed attosccond double-slit interference 
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Initially it was thought there would be a direct corre- 
spondence between the HHG and ATI spectra (see [HI 
and references therein) and attempts were made to ex- 
press the harmonic yield as a sum over ATI channels 
plus recombination [hil. [l5j|. However, no direct link be- 
tween the intensities of individual HHG and ATI peaks 
could be drawn as in general it is not possible to disen- 
tangle the contributions from the different trajectories. 
Two trajectories producing the same harmonic frequency 
will generally lead to different ATI energies. Here we 
show that, for extremely short laser pulses with suitable 
carrier-envelope phase this link turns out to be possible 
because only a very limited number of trajectories con- 
tributes. Taking also advantage of the exponential de- 
pendence of the ionization rate on the field strength, we 
present strong one-to-one links from HHG frequencies to 
ATI momenta, based on shared birth times of the HHG 
and ATI trajectories. We show how to use the relation 
between HHG and ATI to improve molecular imaging 
techniques such as orbital tomography. 

If there is only one trajectory contributing to each ATI 
momentum and HHG frequency, and if both trajectories 
are born at the same time, the ionization steps are iden- 
tical and there is a one-to-one mapping from HHG fre- 
quency u) to ATI momentum p^(uj). We assume here 
that the ATI electron is emitted along the laser polar- 
ization axis, i.e. there is no rescattering ATI. Then, 
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the HHG intensity S(w) = |a(w)|2 and ATI intensity 
A(p) = \M(p)]2 are related by by (atomic units are used 
throughout) 0, 0, El 



state as ipo(v) and using r s = t s — t' s , tq = to — t' , 



a(w) = a(oj) V r ec(w), 
A{ p ^(lj)) =C{u) |a(w)|2. 



(la) 
(lb) 



Here, the complex amplitude ol{lS) is the Fourier trans- 
formed dipole acceleration and a(w) describes the con- 
tinuum wave packet for HHG. The velocity- form recom- 
bination matrix element for the HHG process is denoted 
v rec (w). The factor C(uj) relates HHG and ATI and in- 
cludes the effect of electron motion after the return time 
on the momentum distribution. Below, we confirm that 
C{uj) only depends on the laser field and is independent 
of the atom or molecule. This is in contrast to the quan- 
tity a(w), which is species-dependent. Thus, if the mo- 
mentum distributions A(p) are known for two different 
systems, the ratio of their factors a(oj) can be obtained 
from Eq. (fib]) . 

Before demonstrating the improved molecular imag- 
ing scheme, we find suitable laser pulses for which the 
one-to-one mapping holds. To this end, we express the 
HHG and ATI yields using the strong-field approxima- 
tion and expand around classical trajectories. For HHG 
we employ the Lewenstein model [11] . In this model, 
the saddle-point integration over momentum gives the 
saddle-point momentum k s (t,t') = - f* A(t")dt" /(*-*') 
with A(t) = — f_ E(t")dt" such that an electron born 
at time t 1 returns to its initial position at recombina- 
tion time t. for a linearly polarized laser field E(t). In 
contrast to 0, E3] i we perform both remaining integra- 
tions over t' and t using the saddle-point method. The 
resulting spectrum is generated by trajectories with com- 
plex saddle-point times t' s and t s , starting with imaginary 
initial momentum Vi = i-J2I v and returning with mo- 
mentum v r = ±y2(w — I p ) where J p is the ionization 
potential. We employ a very short pulse such that all 
return momenta v r have the same sign (chosen negative 
here) 18|. The classical times t' , to are defined by set- 
ting Vi = 0. We expand the times t' s , t s around t' Q , t' Q to 
second order in the Keldysh parameter 7, 
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WU%J]' 2 ~ E(t' )(v r + (t -t' )E(t )): 

4C/ P ( E{t ) E'(t' ) 

(E(t' a ))2 \v r + (t -t' )E(t ) E(t' ) 



(2) 

(3) 
(4) 



Expanding to fourth order for the action S(t,t') = 
\ $1 dt" [k s (t, t') + A{t")}2 +I p (t- t'), the resulting ex- 
pression in a d-dimensional world is, denoting the bound 
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+ 0(73), (5) 



with f(t , t' ) 
element 
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r +T E(t a ) 1 E(t' ) 
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The ionization matrix 
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exhibits a pole at the saddle-point momentum for 
Coulombic potentials. Supported by the fact that HHG 
can be modeled succesfully using Gaussian bound states 
that do not exhibit the pole [ll|, we replace the integral 
in Eq. (|6]) by an arbitrary constant. 

Similarly, we expand the ATI amplitude M(p) as given 
by Milosevic et al. (H-Eil around classical birth times. 
For detailed derivations of Eqs. ([2])-([5]) and the analogous 
ATI expression, see [22j |. 

We calculate the uniqueness of a trajectory in deter- 
mining harmonic u by dividing the absolute value of the 
corresponding term in Eq. ([5]) by the total sum. Simi- 
larly, for every harmonic trajectory we also calculate the 
uniqueness of the ATI trajectory born at the same time 
in determining its associated ATI momentum. The maxi- 
mum attainable product of these two factors is labeled 
Q{uj). For the maximum possible value Q = 1, there is 
a perfect correspondence between a harmonic frequency 
and ATI momentum through their shared birth time. In 
Fig.QTa) we show Q(uj) for different two-cycle sin2-laser 
pulses with intensity I = 2x 10 14 W/cm 2 and wavelength 
A = 2000 nm shining on a ID system with the ionization 
potential set to I p = 30.2 eV. The pulses are character- 
ized by the carrier-envelope phase </>cep, i.e., the phase 
between the envelope and the carrier wave of the pulse. 
We only consider ATI momenta whose amplitudes are 
greater than 1 a.u. 

A two-cycle pulse with 0cep = 1.2&7T gives rise to a 
good link between HHG and ATI over a broad harmonic 
range. This result is nearly independent of the dimen- 
sionality d (not shown) . In Fig. QJb) we plot the time- 
dependent electric field and vector potential of this pulse, 
and we indicate the birth and recombination times of 
the dominant trajectories. Because \E(i)\ is much higher 
during the birth time of harmonic orders ~ 100-200 
(t ~ 0.92T) than around t ~ 1.35T— the only other time 
where ATI electrons with the same final momentum are 
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FIG. 1. (Color online) (a) Q(tv) for two-cycle sin2-pulses with 
(f>CEP = 7r (red dashed line), 0cep = 1.25n (black solid line), 
4>cep ~ 1-57T (blue dotted line) and i^cep = l-75n (green dot- 
dashed line); (b) E(t) (x40, black solid line) and A(t) (red 
dashed line) for 0cep = 1.25n. Also indicated are the birth 
(violet plusses) and recombination (blue crosses) times of the 
dominant trajectories. 



born — the link between HHG and ATI arises. Smartly 
selecting experimental phase-matching conditions might 
allow somewhat longer pulses to provide useful links be- 
tween HHG and ATI. However, for slightly longer pulses 
the HHG frequencies are linked to very low ATI momenta 
which may require including Coulomb corrections for the 
classical trajectories. 

We verify the link between HHG and ATI by numeri- 
cal solution of the time-dependent Schrodinger equation 
(TDSE) for ID H+ with varying internuclear distance R. 
We use the softcore potential 



V(x) = 



-Z 



^{x- f)2 + a2 ^{x + f )2 + a2' 



(7) 



where the softness parameter a2 is adjusted such that 
30.2 eV. The TDSE_is_solved on a grid using 

and the bound states 



the split-operator method [23, 
are found by imaginary-time propagation [2fJ . The grid 
length is 24027 a.u. and it contains 143360 grid points. 
After the end of the laser pulse, the wave function is prop- 
agated for two more cycles. The HHG spectrum is cal- 
culated from a windowed Fourier transform of the dipolc 
acceleration and the ATI spectrum is obtained from the 
momentum-space representation of the wave function af- 
ter removing the bound states by windowing out the in- 
ner 40 a.u. in position space. For the ground state with 
Z = 0.731, the HHG spectra and ATI momentum dis- 
tributions for three different internuclear distances are 
shown in Fig. [2j Here we employ the pulse of Fig. Gib). 

Similar to the quantitative rescattering theory 0, [lfl , 
we calculate the ID recombination matrix elements 
free = ('/^(aOMV's^)) using field-free scattering states 
■0s . We obtain numerically exact ij} s by integrating the 
static Schrodinger equation using the Numerov method 
(see e.g. [Ifjj]) on a grid with a total length of 4000 a.u. 
and 320000 grid points. For an electron approaching from 
positive x we set the wave function equal to exp(— ik(uj)x) 
for the two lowest grid points, where k(ui) = \j2(y> — I p ). 
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FIG. 2. (Color online) (a) HHG spectra for ID H^ at 
R = 2.00 a.u. (black solid lines), R = 2.03 a.u. (red dashed 
lines) and 7? — 2.06 a.u. (blue dotted lines), (b) ATI momen- 
tum distributions. The insets show smoothed details of the 
spectra. 
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FIG. 3. (Color online) Vertically aligned | a(c<j) (lines with 
circles) and \J A{p( A *> (to)) (lines with squares) as a function of 
R. We consider harmonic orders 100 (black solid lines), 150 
(red dashed lines), and 200 (blue dotted lines) and their as- 
sociated ATI momenta p (A) = -1.4922 a.u., p (A) = -1.8301 
a.u., and p' A ' = —2.2315 a.u., respectively. The curves were 
shifted vertically for clarity. Also shown are harmonic 150 and 
ATI momentum p' A ' = —1.8301 a.u. from a six-cycle trape- 
zoidal pulse with 1-cycle ramps (brown dot-dashed lines), for 
comparison. 



After integrating upwards, we require that for large pos- 
itive x the wave function is given by 



cil) s {x) = e - ik ^ x + R c e lk(ul)x , 



(8) 



where R c is the reflection coefficient. This leads to a 
normalization constant 



2e- ik ^ x /(ip a (x)+iiP' B (x)/k). 



(9) 



In Fig. [3] we demonstrate the link between HHG and 
ATI as a function of R using the other parameters of 
Fig. [2] For three different harmonic orders we plot |a(w)| 
and \J A(p( A } (w)) (both normalized) using the links be- 
tween HHG and ATI obtained from Eqs. ([5]) and the 
corresponding ATI expression. The remarkable overlap 
between HHG and ATI that only breaks down for the 
largest R confirms the strong link between HHG and 
ATI indicated by the black solid line in Fig. Ula). The 
ratios of normalization constants correspond to \J C(co) 
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(see Eq. ([lb])). The oscillation in R of both the HHG and 
ATI amplitudes in Fig. [3] is caused by excitation to the 
first excited state before ionization. Between two maxima 
the first excited state drops in energy by exactly ui. We 
have verified that the oscillation period is inversely pro- 
portional to the laser wavelength prior to the moment of 
ionization. Both the fact that the HHG and ATI curves 
for a six-cycle pulse are less similar to each other and 
the fact that they exhibit wild behavior are related to 
multiple trajectories contributing to the yields. 

The link between HHG and ATI gives the experimen- 
talist access to the ratio of the instanteneous ioniza- 
tion rates of different molecules during the high-harmonic 
generation process, and as such is a useful tool in study- 
ing HHG and molecular imaging. In particular, the esti- 
mate for the continuum wave packet needed for the to- 
mographic reconstruction of molecular orbitals Q can be 
improved using 

ag(u) = oW (W) ^Ag (p(A) ( W ))/A(») (p(A) (a,)) . (10) 

Here is the orientation of the molecule in the laser 
field and with the superscript '(a)' we denote quanti- 
ties belonging to the reference atom in the tomographic 
procedure. Demonstrating numerically the possibility of 
combined HHG-ATI molecular imaging, we retrieve the 
field-free matrix elements of the first excited state of ID 
using a reference atom. We employ 



Wroc(w) — a(cj)/a(u)) 



(11) 



with cither a(ui) = {u)F\/P^> (HHG imaging), 
where Pi is the total ionization probability, or with 
Eq. (JlOj) (HHG-ATI imaging). Here we use for ID H^ 
the parameters R = 2 and Z = 1.3 a.u. and the same 
laser pulse as for Fig. [2j As the reference atom we use 
ID softcore models with different nuclear charges Z^ a \ 
For all systems the softcore parameter a2 is adjusted such 
that I p = 30.2 cV. The results of the simulation can be 
found in Fig. S) The figure shows that when the total nu- 
clear charge is identical for the molecule and the reference 
atom (Fig. SJb)), the molecular matrix element can be 
accurately retrieved using only HHG from the molecule 
and atom [22I, thereby also demonstrating the accurate- 
ness of our field-free matrix elements 0. However, when 
the total nuclear charge does not match (Figs. Hta) and 
IDJc)) the propagation step becomes different for the atom 
and molecule and errors arise in the retrieval of the ma- 
trix elements. These errors largely disappear by incor- 
porating also ATI electrons in the retrieval procedure, 
demonstrating the potential of Eq. (TTU)) for orbital tomo- 
graphy and molecular imaging in general. The shallowing 
of the retrieved matrix elements comes from diffusing the 
HHG and ATI spectra with Gaussians with 1/e- widths 
of Aw = 6cjl a.u. and Ap = 2-\/0.2wl a.u., respectively. 

In summary, we used the saddle-point approximation 
and expansions in 7 to evaluate strong-field expressions 



for HHG and ATI in terms of sums over classical trajec- 
tories. Using these expressions we have shown that for 
extremely short laser pulses and long laser wavelengths 
there exists strong links between individual frequencies 
and momenta of HHG and ATI. We demonstrated these 
links and their potential for molecular imaging using ID 
model calculations. Future molecular imaging experi- 
ments will benefit from this effect. 
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